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Abstract 

O , We introduce new discretizations of the action for static quarks. They achieve an 

exponential improvement (compared to the Eichten-Hill regularization) on the signal to 
noise ratio in static-light correlation functions. This is explicitly checked in a quenched 
^ ■ simulation and it is understood quantitatively in terms of the self energy of a static quark 

and the lattice heavy quark potential at zero distance. We perform a set of scaling tests 
in the Schrodinger functional and find scaling violations in the 0(a) improved theory 
to be rather small - for one observable significantly smaller than with the Eichten-Hill 
regularization. In addition we compute the improvement coefficients of the static light 
axial current up to 0(5(q) corrections and the corresponding renormalization constants 
non-perturbatively. The regularization dependent part of the renormalization of the 
b-quark mass in static approximation is also determined. 
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1 Introduction 

Lattice QCD is playing an important role in the interpretation of B-physics experi- 
ments [1]. However, as new physics is hiding behind the standard model, the required 
precision for lattice computations of B-meson transition amplitudes is becoming more 
and more demanding. The development of new algorithms [2-5] promises that light dy- 
namical fermions at small lattice spacings, a, can be reached in the near future, which 
will reduce one of the most important systematic errors. 

Still, the treatment of b-quarks on the lattice remains another difficult part of these 
computations, since it appears unlikely that lattice spacings small enough to satisfy 
a < 1/my, will soon come into reach. We refer the reader to reviews for an explanation 
of the various approaches to this problem [6-8]. 

A theoretically clean solution is provided by HQET. This effective theory starts 
from the static approximation describing the asymptotics as mb — > oo. Corrections 
0(l/mb) have to be computed by a 1/mb expansion, where the higher dimensional 
interaction terms in the effective Lagrangian are treated as insertions into static corre- 
lation functions. An attractive feature of this approach is that the continuum limit exists 
and results are independent of the regularization. The theory requires non-perturbative 
renormalization, but a concrete way to carry this out has been proposed and tested for 
a simple case [9] . 

Furthermore also just the leading order term (static) is of considerable interest, 
since it provides a limit of the theory, which in many cases is not expected to be far 
from results at the physical point. Other methods to treat heavy quarks can thus be 
tested by checking whether they smoothly approach the static limit 1/mb ~^ 0. In 
fact, results can also be obtained from interpolations in 1/mb between points below the 
b-quark mass and the static limit, see [10] for a precise demonstration. 

As we will discuss in detail in section 12.11 the noise-to-signal ratio of static-light 
correlation functions grows exponentially in Euclidean time, ~ exp(^XQ) and the rate 
// ~ 1/a diverges as one approaches the continuum limit. Good statistical precision is 
difficult to reach for xq > Ifm [11,12]. In the past, sophisticated wavefunction tech- 
niques have been used to obtain ground state properties [13-15]. The results of these 
attempts, obtained with the standard Eichten-Hill regularization [16], were not com- 
pletely satisfactory in all cases [17]. In [18] we used the fact that fi is not universal. By 
changing the regularization, we were able to obtain results at large xq , where the ground 
state can be isolated with good statistical and systematic precision. Here we will discuss 
these alternative actions in more details. Particular emphasis is put on the size of dis- 
cretization errors, since their smallness should always be the first criterion for choosing 
an action. Let us mention right away that we will consider only actions, which share the 
symmetries, eqs. ()2.18I2.19() . with the Eichten-Hill action, since these guarantee that no 
0(a) terms are needed in the static action to have 0(a) improvement. For all actions 
considered, we determine the coefficients of the improvement terms needed for the axial 
current with 1-loop precision. In particular we correct for a mistake in the perturbative 
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computation of the improvement coefficients in [18] , see the Erratum. Furthermore we 
determine the action-dependent parts of the renormahzation of the static axial current 
and the b-quark mass. 

Before coming to the description of the static actions, we hst some preliminaries. As 
a probe of the theory, we will consider correlation functions defined by the Schrodinger 
functional (SF) [19,20]. For the introduction of static quarks in the SF as well as 
any unexplained notation we refer to [21]. We adopted the 0(a)-improved Wilson 
regularization [22,23] with non-perturbatively determined value of csw [24] for the light 
quarks. The gauge sector is discretized through the Wilson plaquette action. 

Below we will report also on results of Monte Carlo computations. Our simulation 
parameters are summarized in appendix O All the results have been obtained in the 
quenched approximation, a part of them already appeared in [18] and [25]. 



2 Static actions 

Static quarks are introduced on the lattice through fermionic fields Vh and ip^^, which 
live on the lattice sites and satisfy the projection properties 

1 + Tn 

P+V^h = ^h, Vh^+=Vh, P+ = ^^. (2.1) 

The action Sh for static quarks has been derived by Eichten and Hill in [16]. Here we 
adopt a notation, which allows us to provide a generalization of that action and write 
it in the form 



= E V'h(^)(^o'^ + Smy,)M^) , (2.2) 

1 + a omw ^ 

with the covariant derivative 



D^il^l^ix) = i il^hix) - W^{x - aO, 0)V'h(2; - aO) , (2.3) 

where W{x,0) is a gauge parallel transporter with the gauge transformation properties 
of the link U{x,0). The Eichten-Hill action is given by setting W{x,0) = U{x,0). The 
quark propagator {x^y) satisfying {D^ + 6m^^)G'^ {x,y) = 6{x — y)P-{-, reads 

G^{x,y) = 9ixo - yo) <5(x - y) (1 + a <5mw)-("""^°^/'^ r'^iy,x)^ P+ , (2.4) 
V'^ix,x) = l, (2.5) 
P^(x, x + Rfi) = W{x, n)W{x + afi,fM)... Wix + a{R - l)fi, fi) for R > 0, (2.6) 

where 5my/ cancels the divergence in the self-energy of the static quark (see below). 
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2.1 Noise to signal ratio 

In order to simplify the following discussion we start from the action in eq. (|2.2)) but 
set (5mw = 0. This is possible since the dependence of correlation functions on 5mw is 
known exactly from the form of the static quark propagator in eqs. 1)2. 412. 6() and can be 
restored easily. The noise to signal ratio discussed here is independent of Sm^j^. As an 
example we consider a correlation function of heavy- light fields such as f^^^ixo), which 
describes the propagation of a static-light pseudoscalar meson. Its definition 



/r(^o) 



a 

y 



5](Ar(x)Ch(y)75Ci(z)) 



(2.7) 



involves the boundary fields and Ci (see [21]) as well as the time component of the 
axial current 



{21 



where the subscript '1' indicates the light quark fields. The correlation function is 
represented pictorially in figure ^ From the quantum mechanical representation one 




Xn=0 



Xn=0 



Figure 1: Schematic representation of the correlation functions /A^*(a;o) and in the SF. The 

double lines indicate the static quark propagators. 



expects the large xq asymptotic behavior 



(2.9) 



where -Egtat is the binding energy of the static-light system. It depends on the choice 
for W and diverges approximately linearly in the lattice spacing. 



Est^t ~ ^seif + 0(a°) ~ ie(i) 50 + • • • • 



(2.10) 



This divergence is to be canceled by 5mw = — -E'seif + 0(o*^). The finite part of (5mw 
is of course scheme dependent, but in perturbation theory one can always perform the 
double expansion 



^self = E^l},gl + 0{gl) 



E. 



(1) 

self 



(2.11) 
(2.12) 



For example, if we define -Egoif from a Schrodinger functional correlation function (see 
eq. ()A.13|) ). the a-expansion appears as an expansion in terms of a/L. The coefficient 
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e^^^ does then depend on the action, but not on the correlation function used to define 
it. 

We now want to discuss the noise to signal ratio of the Monte Carlo estimate of f^^^. 
The variance of a generic observable O in a QCD simulation is given by (0^)g ~ (^)g' 
where the expectation value (...)g is an average over the gauge fields including the 
weight obtained after integrating out the fermion fields (the fermion determinant). O is 
derived from O by performing the corresponding Wick contractions on each gauge field 
background. In the Monte Carlo one takes advantage of translation invariance on the 
r.h.s. of eq. ^ and estimates P^^^ = (O) with O = Ey,z,x ^r*(^)Ch(y)75Cl(z)- 

Its variance, aAi^o), can be rewritten as 

fTA(xo) = ^ Yl ([^r*]hl(a;0,x)Ch(y)75Cl(z)x 
x,y,z,x',y',z' 

{[<^\,i'(xo,x')Ch'(y')75Ci'(z')}^) - [/r(^o)]^ (2.13) 

where h' and 1' denote copies of h and 1 differing only by their flavor, which are introduced 
to be able to write the variance in the form of a standard expectation value. Saturating 
crA(a;o) with intermediate states shows that its large xq asymptotics is 

aA(xo) oc J2 C(x,x')e-^«^"'('^'''') , (2.14) 

x,x' 

with Eii'{x, x') the energy of a state with a static quark-antiquark pair at positions x, x' 
and a light quark-antiquark pair with flavors 1,1'. At large xq the sum is dominated by 
the term with the smallest energy E. This is expected to be given by x = x', where 
the color charges of the static quark and anti-quark compensate. The light quark- 
antiquark pair should then feel little of the static quarks and its lowest energy be given 
approximately by m^. This argument suggests 

Minx,x' En'ix,x) = V{0) + . (2.15) 

Here V{0) is the lattice static quark potential at zero distance, which for example is 
computable via the large xq asymptotics 

(?h(^)75V'h'(^) Ch'(x)75Ch(x) ) oc e-^«^W . (2.16) 
The noise to signal ratio i?NS for /^*^*(xo) should then approach 

iiNSOce[^---(™-+^(°))/2]x.o _ (2.17) 

For the Eichten-Hill action (or more generally whenever W{x,0) is unitary), the lattice 
potential vanishes at zero distance and one obtains the formula given earlier by Lep- 
age [6] . We infer from eq. 1)2. 17() that the linear divergence in -Estat is responsible for the 
exponential growth of the error with the Euclidean time xq, which has been observed in 
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many numerical investigations of correlation functions such as fj^^^ [12-15]. It is then 
clear that the problem becomes more severe as the continuum limit is approached. In 
particular, for the Eichten-Hill action the coefficient e^^^ was found to be rather large [26] 
rendering precise computations hopeless. Finally we remark again that the form of the 
static propagator shows that i?NS is (exactly) independent of 5mw- In the final formula 
eq. 1)2. 17(1 this is realized since Sstat and V{0)/2 are shifted by the same amount, 

aEste^tlsmw = «^statl5^^=o + ln(l + a(5mw) 
«^(0)l5„w = «V^(0)l5™^=o + 21n(l + a<5mw). 



2.2 Statistically improved discretizations 

A possible way to improve on the problem discussed in the previous section is to employ 
actions inspired by variance reduction methods. Such methods led for example to 
the introduction of the one-link integral (or multihit) [27] in the pure gauge theory. This 
provides unbiased estimators for quantities such as Polyakov loop correlation functions, 
significantly reducing at the same time their variance. Here, a similar change will not 
lead to an unbiased estimator, but rather has to be considered a change of the static 
action. Alternatively (or simultaneously) one can try to enhance the signal by choosing 
the parallel transporter W such that £^stat ~ -^seif is comparatively small (assuming 
V{0) to be numerically less relevant). In addition to this, we want to preserve on the 
lattice the following symmetries of the static theory 

i) Heavy quark spin symmetry: 

tph — > Viph , iph — > i'h^'^ > with V = ex.p{-i(l)ieijkcrjk) , (2.18) 

ii) Local conservation of heavy quark flavor number: 

^ e^^W V'h, V^h ?he-^''W. (2.19) 



Together with gauge invariance, parity and cubic symmetry, this is enough to guarantee 
that the universality class and the 0(a) improvement are unchanged with respect to 
the Eichten-Hill action, as has been discussed in [21]. Furthermore we want to keep the 
action as local as possible. We therefore exclude constructions of W which involve fields 
at a distance two lattice spacings or more away from the link U{x,0). Taking these 
considerations into account, we propose the following regularized actions 



QS 



c-HYP 



W^{x,0) = 
W'{x,0) = 



Vix,0) , 
V{x,0) 

^HYP(X,0) 



^+ (ltTV\x,0)V{x,0) 
5 \ 3 



1/2- 



(2.20) 
(2.21) 
(2.22) 
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where V{x,0) is the average of the six staples around the Unk U{x,0) 

1 ^ 

^(^,0) = -Y,[u{x,j)U{x + aj,0)U\x + a6,j) 
i=i 

+U^{x-aj,j)U{x-aj,0)U{x + a6-aj,j)\ , (2.23) 

and Vhyp is the HYP hnk [28]. In the latter case three coefficients {ai,a2,as) = a 
need to be specified in order to define the combination of differently smeared links 
in the construction of the HYP link. In the following we will only discuss the choices 
a = (0.75, 0.6, 0.3), motivated in [28] and a = (1.0, 1.0, 0.5) obtained by an approximate 
minimization of Rj^s- They define 5^^^^ and S"^^^^ respectively. The construction 
of HYP links involves projecting 3x3 complex matrices onto SU(3). As we expect 
deviations from SU(3) to be small, we prefer approximating the projection by an analytic 
function defined by the steps 

W / ^JtI{WW^)/3 , (2.24) 

followed by 4 iterations of 

W ^X{l-^Im (det X)) , with X = W - ^W^w) . (2.25) 



To be precise, this function is to be taken as part of the definition of VhyP) when our 
numerical results for improvement coefficients and renormalization factors are used in 
subsequent computations. 

For the actions S^^^^ and S^^^'^ the value of e*^^^ has been computed in [29]. In 
the cases W ^ SU(3) it is more convenient to look at itself — ^(0)/2 not only because 
this is the relevant quantity for the noise to signal ratio, but also because at 1-loop order 
the deviations of the links W from unitarity cancel in the combination. We indicate by 

the 1-loop coefficient in the perturbative expansion 

E,,u - V{0)/2 ~ Qr(i) + 0(aO)) gl + 0{g^) . (2.26) 

Of course r^^^ = e^^-* for S^^^ and S^^. In appendix ^ we outline a computation of 
for and explain why it is the same for (at one-loop order). Some results for 
r^^^ are collected in table ^ From the table we see that perturbation theory already 
suggests S^^^^ as the most favorable choice concerning signal enhancement. 

The discretizations and are inspired by noise reduction methods, namely 
APE smearing and one-link integral respectively. APE smearing was introduced in [30] 
where it was shown to suppress fluctuations in gauge invariant quantities in the pure 
gauge theory. As shown by this can also be interpreted as a reduction of the 
static self energy. Finally is an approximation of the SU(3) one-link integral. For 
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aV{0) 




0.16845(2) 


0.68(9) 


0.0 




0.05737(4) 


0.85(2) 


0.671(2) 


OS 


0.05737(4) 


0.76(2) 


0.496(2) 


^HYPl 


0.04844(1) 


0.44(2) 


0.0 


^HYP2 


0.03523(1) 


0.41(1) 


0.00001(1) 



Table 1: One loop coefficients r^^\ defined in eq. J2.26'|) and non-perturbative values for aiSstat and 
aV{0) at 13^6. 



completeness we describe briefly how we arrived at in appendix^ The reader may, 
however, take it just as another ansatz satisfying our criteria explained above. 

The effectiveness of these regularizations in reducing the noise to signal ratio in 
static-light correlation functions has been checked non-perturbatively by a simulation 
on a 16^ X 32 lattice at P = 6/(7q = 6, where the lattice spacing is a ~ 0.1 fm (the 
hopping parameter k has been set to the strange quark mass value, k. = 0.133929, as 
in [18]). Figure 121 shows the results for Rf^s obtained using an ensemble of almost 5000 
configurations. 



R 
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10 20 



30 



on a 16'' X 32 lattice at (3 — 6. Filled circles refer to S™, empty circles to (Sl^ gives similar results), 
empty (filled) triangles to S^^^^ (S^^^'^)- 



The largest improvement is again given by the action 5^^^^, but we see from the 
figure that all the proposed discretizations produce a clear signal (for the considered 
statistics) at least up to a time separation of roughly 2 fm. The dotted lines in figure [21 
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represent the predictions from eq. ()2.17|) . where for aE'stat; aV{0) and aruT^ we insert 
the estimates from the data. We smnmarize the obtained values of a-Egtat and 0^(0) in 
table E The formula in eq. (|2.17j) turns out to be always quite accurate in describing 
the results. Furthermore, assuming -Estat — ^(0)/2 to be dominated by its divergent 
term and approximating it by the leading perturbative estimate yields the correct order 
for -Rns- It is thus to be expected that the ordering of i^NS observed in figure HI will 
be preserved when going to smaller lattice spacings, where it is even more important 
to keep i?NS at a reasonable level if one wants to reach distances xq around 1 — 2 fm. 
Indeed, numerical experience supports this expectation [18]. 

3 The size of discretization errors 

In addition to the reduction of statistical errors, the size of scaling violations has to 
be taken into account when one chooses between different actions. For this reason we 
designed a number of scaling tests by which we could study the approach of a set of 
quantities defined in the SF to their continuum limit values. Before going into the 
details of this study we first need to discuss the computation of the 0(a) improvement 
coefficients of the static axial current for the different actions introduced. 

3.1 Computation of the improvement coefficients cj*^* and 6^*^* 

The 0(a) improvement programme [31, 32] has been carried out for the EH action 
in [21]. In [18] the discussion has been extended to the actions we are considering here. 
As mentioned in the previous section, arguments based on the symmetries preserved 
on the lattice allow to show that the improvement pattern is the same in all cases. In 
particular the actions are improved once the light sector has been improved; no new 
improvement terms are needed for the static part of the actions. 

Concerning the static-light axial density ^q*^*(x), the improved version reads 

«^*)o = A«o*'*(^)+acr'^<'*(^) , ^K"' (x) =Mxhjl5 ^'^^*' V'h(x) , (3.1) 

where the improvement coefficient cj*^* has been introduced. In a mass independent 
renormalization scheme the renormalized density can be written 

{At% = Zt\go, a/x)(l + bt'am^){Af-% , (3.2) 

with mq the bare subtracted light quark mass (cf. eq. HC.1|) '). Zpf^{gQ,afi) the scale 
dependent renormalization factor of the axial current and 6^^* a second improvement 
coefficient. The values of (f^^ and 6^*^*, as well as Z^*^*, depend on the choice for the 
static action. For the EH action the improvement coefficients cj'^* and 6^^* have been 
determined at 1-loop order of perturbation theory in [21,33,34], while the renormaliza- 
tion constant Z^^^^ has been computed non-perturbatively in the SF scheme in [35]. In 
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the remainder of this section we describe our computation of the improvement coeffi- 
cients cj'^* and 6^^* for the actions in eqs. (|2.2UI - 0.22() . and we present a set of scahng 
studies. We will come back to the renormalization constant Z'j^^^ in section 4, where 
we also discuss the improvement the new actions can bring in the computation of the 
b-quark mass following the strategy in [9,36]. 

3.1.1 Improvement conditions and results for cj'^* and 6^^^* 

We want to compute the 1-loop coefficients c^'^*''-^^ and 5^^*'^^^ of the improvement 
constants c^*^* and 6^*^*. To this end we have adopted a mixed strategy. For and Sf^ 
the computation has been carried out analytically^, while for S^^^^ and 5^^^^ we have 
used Monte Carlo simulations to numerically estimate an effective 1-loop coefficient for 
the coupling range relevant here. In both cases we have exploited the same improvement 
conditions. We collect some details of the analytic computation in appendix \K\ 

We introduce the correlation ff^f defined as ff^'^^ in eq. (|2.7|1 with the current ^q**^* 
replaced by the improved current in eq. (|3.1j) . Assuming the knowledge of cj^* for a 
discretized action S^, we have enforced the condition^ 



r^fiT/2,Sl) ^ ftf{T/2,Sl) 



ftf{T/2,Si) ftf{T/2,Sl] 



e = 0, 9' = 1 , mq = , (3.3) 



and solved the implicit relation for c^*** of S^. Expanding eq. ()3.3|) in powers of the 
coupling 5q and setting = S™, for which c^^*'^^^ is known from [34], we obtain 

4tat^^^at,(i)^2^Q^^4)^ c^^-^'W = 0M72{4) , for 5^ and 5^ (3.4) 

The ratios in eq. (|3.3j) have then been evaluated, with = S^, on the configurations 
generated in runs I to IV of table El in order to obtain cj*** for = S^^^^, S^^^"^. 
There the size L/a and the /3 values have been chosen such that L is kept fixed to 1.436ro 
with ro from [37,38]. Our results are shown in figure 01 The result for provides a 
test of the numerical procedure since it is known analytically (eq. (|3.4|) ). We see that 
by ascribing to the numerical result an error, which covers the spread of the points we 
find agreement with eq. 1)3.4(1 . This suggests that here higher orders in Qq contribute 
little to the cutoff effects of the ratios appearing in eq. (|3.3|) . We estimate in the same 
way the errors for cj'^* of 5^^^^ and 5^^^^, quoting the result 

cT' « 0.039(4) , for ^J^^Pi (3 5) 

cT' ^ 0.220(14) 5^ , for S^^^^ . (3.6) 



^At this order in perturbation theory the results for the improvement coefficients are the same for 
the two regularizations, see appendix IXI 

^Where necessary we explicitly indicate the dependence of the correlation functions and the improve- 
ment coefficients on the discretization of the static action. 
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0.2 



w o 
ClO 



0.1 



-0.1 



QHYP2 



QHYPl 




0.92 



0.94 



0.96 



0.98 



Figure 3: Numerical results for the improvement coetEcient c^^*. For the result of [34] is plotted, 
for the band represents the result eq. 13.411 . while otherwise it refers to the numerical estimates 

eqs. nun:^ . 



Here the statistical errors are not the relevant ones. Rather our dominant error is due to 
the assumption (tested to some extent as just explained) that the 1-loop value for cj'^* 
for action = is accurate for our values of g^. Apart from this, our determinations 
do provide a non-perturbative computation of cj*^* for the other actions. Thus the errors 
quoted in eqs. 1)3. 513. 6() include a reasonable (but not necessarily a safe) estimate of the 
perturbative uncertainty. 

Note that for S^^^^ the improvement coefficient is considerably larger than for the 
other actions. This means that before improvement the linear a-effects are larger in this 
case and one might expect that this will be the case also for higher order a-effects. We 
will see below that this is however not born out of our non-perturbative results. 

Sensitivity to the improvement coefficient b^^^^ is achieved by exploiting the quark 
mass dependence of the correlation function f^f'ixo). Again, supposing 6^*^* be known 
for the action S^, we consider the improvement condition 



{l + abr\Sl)m^,)fff{T/2,Sl) 



{l + abT{Si)m^')fff{T/2,S, 



{l + abt\Sl)m^)P^f{T/2,Sl) 



[l + abt\Sl)m^)P^f{T/2,Sl) 



(3.7) 



with 6 = 1, and solve for ftj'^* of S^. 
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In perturbation theory 6^^* has been computed to 1-loop order in [21] for the 
Eichten-Hill regularization. Similarly to what we did for cj'^*, we have then analytically 
computed 6^^* for S^, i.e. we have expanded eq. (|3.7|) in perturbation theory with 
= 5™ and 5*^ = 5^. In this computation we set ruq = and, following [39], mq/ 
such that Lm^i = 0.24, with m-R/ the quark mass renormalized at scale fi = 1/L in the 
minimal subtraction scheme on the lattice. The result reads 



bt' = 1/2 + 6f 5o + 0{g^) , 6f = 0.033(7) , for SI and . (3.8) 



Next we impose eq. p.7|) on the data sets I to IV in table 1^1 Here, both k values and 
6 = 1 are used, where the second k value was determined such that Lm-^i = 0.24 in the 
SF-scheme at scale /x = l/(1.436ro). With = S^, we obtain 



T,stat 
T,stat 



1/2 + 0.078(8) gg, 
1/2 + 0.259(13) 5^ 



for 



for 



(3.9) 
(3.10) 



The errors have been estimated as in the case of c^***. We show the numerical results 
in figure Again the improvement coefficient turns out largest for 5^^^^. 

0.3 



CV 

T 1 

I 



m <; 



0.2 - 



0.1 - 



-0.1 




0.92 0.94 0.96 0.98 

Figure 4: Numerical results for the improvement coeScient 6^^*; presented as in figure]^ 

One comment is in order here. In the condition in eq. (|3.7j) . 0(amq)-terms with 
coefficients bg and bf^^^ [21] are neglected, while terms proportional to b^ (see [40]) drop 
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out. Of course, is 0{gQNf) and therefore irrelevant in a 1-loop computation. While 
bg is 0{gQNi), it enters expectation values only with an additional factor (unless there 
is a non-zero background gauge field). In other words, eq. (|3.7|) is correct up to 0(5q) 
in full QCD, but in the quenched approximation it is valid also non-perturbatively. 



3.2 Scaling tests 

On the same set of configurations used for the determination of cj*^* and 5^*^* (runs 
I-IV), and for k = Kc, we have computed (for each regularization) the ratios 



/iT(r/2) 



Here and in the following f^^^{xo) refers to the 1-loop improved axial current in eq. (|3.1j) . 
The correlation function ff^^^ is schematically represented in the right part of figure^ 
it precisely reads 

/r = E (Ci'(u)75Ch'(v) Ch(y)75Ci(z)) , (3.12) 

u,v,y,z 

with the primed fields living on the boundary at xq = T. The quantity h[d/L) was 
introduced in [35]. It is defined through the boundary to boundary SF correlator 

/i'\^3) = -^ E (Ch'(x)75Ch'(0)Ch(y)75Ch(z)), (3.13) 

note that the sum runs on xi and X2 and therefore yields an X3-dependent correlation 
function. 

All the quantities in eq. (|3.11|) have a finite continuum limit in a fixed, finite, 
volume. The continuum result is universal, but the way it is approached depends on 
the details of the regularization. The scaling of these ratios can then provide an idea 



about the size of discretization effects for the static actions in eqs. (|2.20l - 12.22)1 . Our 
results are shown in figure El The cutoff effects in are somewhat larger with the 
Eichten-Hill action than with any of our alternatives. Still, it is more relevant to note 
that all cutoff effects visible in the figure are very small and linear in a?'. This suggests 
that the 0(a) improvement programme has been implemented in a satisfactory way for 
the lattice spacings considered here - also for Sj^^^^, where the improvement coefficients 
turn out to be not that small. 

The gain in statistical precision brought by the new discretizations of the static 
action can also be seen in figure especially for the quantity /i(l/4), which involves 
two static quarks propagating over the whole temporal extent. 

Another interesting observable is the step scaling function introduced in [36] for 
the renormalization of the b-quark mass and further discussed in the following section. 
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0.015 



(a/L)^ 



Figure 5: Scaling plots for £,a(0, 0.5), ^i{0, 0.5) and /i(l/4). Symbols are as in figure0 some of them 
have been shifted on the horizontal axis for clarity. Again the results for are similar to the ones from 
and haven 't been plotted. 



Here we simply define it as 



S^(n,a/L) = 2L[rstat(2L) -rstat(i)]„= 



(3.14) 



rstat(^) = ^ln[/i7(xo-a)//if(xo + a)] , = 1/2 , T = L , (3.15) 

in terms of the correlation function /^*f* and the Schrodinger functional coupling g{L) 
[19,41,42], which fixes the length scale L. The continuum limit of T,^ is, for each value 
of the coupling u, a universal quantity, independent of the regularization used. The 
step scaling function can hence be used for a further scaling study, in particular since 
discretization errors with the EH action turn out to be clearly visible [9]. Figure IHl shows 
that they are rather linear in (a/L)-^ for all actions considered. Compared to the EH 
action, the alternative ones show smaller lattice artefacts. In fact they are the smallest 
for S^^^"^ which also showed the best noise to signal ratio. We will return to another 
scaling test in section ITTl 
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Figure 6: The step scaling function Em for the different static actions. Symbols are chosen as in 
figure^ Results for the Eichten-Hill action are taken from [9]. 



4 Renormalization of the axial current and quark mass 

In the effective theory the static-hght axial current is not derived from a symmetry 
transformation of the action. The renormalized current is therefore scale dependent. 
This dependence has been studied non-perturbatively, over a wide energy range, in the 
SF scheme (see [35]). The main quantity considered in that study is the step scaling 
function SJ'^*(m, a/L) defined as 



:t\u,a/L)- ^r(50,2L/a) 



Zf-'{go,L/a) 



(4.1) 

g^{L)=u, mq=0 



with 



.Stat , ^ fffiL/2) 



(4.2) 

e=o.5 



here H^'^^ is the tree-level value of S, as in [35]. In eq. ()4.2|) fi is the correlator between 
two light-quark pseudoscalar boundary sources 

/i = -|^ E (Ci(u)75C^(v)C2(y)75Ci(z)), (4.3) 

u,v,y,z 

and we remind the reader that g{L) is the Schrodinger functional coupling. In the first 
part of this section we report on a scaling study of S J'^* . 

In addition, for a chosen reference scale Lj-ef = 1.436ro, we are interested in the 
dependence of the renormalization factor Z^'^^(gQ, Lref/a) on the bare coupling gQ. This 
allows to match bare matrix elements of the axial current to continuum ones (up to 
cutoff effects). The renormalization factor itself clearly depends on the regularization, 
and needs to be recomputed for the actions in eqs. (|2.2()I - [T^^ . 

In the quenched approximation the b-quark mass has been obtained at leading 
order in HQET through the matching of the effective theory to the full one [9,36]. The 
mass of the B-meson is used in the procedure as phenomenological input. Therefore, 
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although the matching is performed in a small volume, the effective theory has then 
to be connected to large volumes. Improving on this part of the computation would 
reduce the error on the result for the b-quark mass. The use of the static actions 
introduced here (in place of the Eichten-Hill action used in [9,36]) has been proven to 
be very efficient in this respect (see [8], where S^^^^ has been considered). In the last 
part of this section we briefly discuss the strategy which has been used to compute 
the b-quark mass, with particular emphasis on the quantities we recompute in the new 
discretizations. 



4.1 Results for the step scaling function and Z'^'^^ 

The step scaling function in eq. 1)4. 1() has been evaluated for u = 3.48 and with lattice 
resolution L/a = 6,8, 12 on the configurations produced in the runs ZI-sZIII of tableEl 
(the same sets were already used in figure El . In figured the results are compared to 
those for the Eichten-Hill action, obtained in [35]. In this case a larger set of lattice 
spacings was used and the data have been extrapolated to the continuum limit. Here 
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Figure 7: The step scaling function E(3.48, a/L). Symbols are as in Egure\^ some of them have been 
shifted on the horizontal axis for clarity. The results for 5*^ are similar to those from and haven't 
been included in the plot. Data for the Eichten-Hill action and the continuum limit value have been 
taken from [35]. 



again we see that the different regularizations yield very similar results already for finite 
lattice spacing, although in principle agreement needs to be found only after having 
taken the continuum limit. 

At the scale L^f = 1.436ro we have also computed the renormalization constant 
^stat ^YiQ /3-range relevant for large volume simulations (again data from runs I-IV 
have been used here). We summarize the results in table |2j For 6 < /? < 6.5 we 
parameterize them in the form (x = /? — 6) 



Zt^\go,L,,,/a) 
Zt\go,Lref/a) 

Zt\go,Lref/a) 



0.7889 - 0.1290 x-h 0.1197 
0.7895 - 0.1291 0.1189 
0.7962 - 0.1368 x + 0.1207 
0.8443 - 0.1740 x + 0.1351 



for 

for SI 
for 5^^' 
for 



(4.4) 
(4.5) 
(4.6) 
(4.7) 
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p 


qA 


OS 


= 

qHYPI 


QHYP2 


6.0219 
6.1628 
6.2885 
6.4956 


0.7873(9) 
0.7722(10) 
0.7633(10) 
0.7565(12) 


0.7880(8) 
0.7727(10) 
0.7640(10) 
0.7569(12) 


0.8015(7) 
0.7845(9) 
0.7746(9) 
0.7653(11) 


0.8823(7) 
0.8579(9) 
0.8409(10) 
0.8227(11) 






„stat 
^A — 


1-loop 

cHYPl 


QHYP2 


6.0219 
6.1628 
6.2885 
6.4956 


0.7862(8) 
0.7708(10) 
0.7619(10) 
0.7543(13) 


0.7868(8) 
0.7714(10) 
0.7624(10) 
0.7547(13) 


0.7933(7) 
0.7770(9) 
0.7669(9) 
0.7580(12) 


0.8406(7) 
0.8196(8) 
0.8054(8) 
0.7913(11) 



Table 2: Results for ZX""* at the scale L^cf = 1.436ro. 



when c^*^* is set to its 1-loop value, and 



^A 

'7Stat 
^A 



{go,Lrcf/a) 
igo,Lrcf/a) 
Zt\9o,Lref/a) 
Zt\go,L,,,/a) 



0.7900 - 0.1286 x + 0.1231 
0.7907 - 0.1287 x + 0.1224 
0.8044 - 0.1404 x + 0.1244 x^ 
0.8866 - 0.1995 x + 0.1423 x^ 



for St 



for Si 



HYPl 



for 5l 
for 5^^' 



(4.8) 
(4.9) 
(4.10) 
(4.11) 



for (f^^ = 0. The formulae reproduce the numbers in table |21 within their errors. They 
may be assigned an uncertainty of about 2%o. 

Given now a bare matrix element $bare of the static-light axial current, the corre- 
sponding matrix element $(//) renormalized in the SF scheme at the scale //~^ = 1.436ro 
can be written 



$(^) = Zi*^*(c/o,ircf/a) ^bare(50) , ^^ ^ = Lrel , 



(4.12) 



with Z^'^^{gQ, L^Qf/a) from eos. 1)4.4114. 11|) depending on the regularization used to com- 
pute the bare matrix element. Finally the Renormalization Group Invariant (RGI) 
matrix element $rgi reads 

"^^^^ ci>(^) , (4.13) 



$RGI 



$(/.) 



where we have factorized out the universal ratio $rgi/'^(/^)) which has been computed 
in the continuum limit in [35]. There one finds the result 



$RGI 



1.088(10) , at = (1.436ro)~^ 



(4.14) 
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which, as discussed, applies also to the discretizations introduced here. Moreover, as 
the error on the ratio in eq. H4.14() refers to a continuum result, it propagates into $rgi 
only once has been extrapolated to the continuum limit. 

4.2 Renormalization of the quark mass 

In [9] a strategy for the non-perturbative computation of the b-quark mass in the static 
approximation has been introduced and discussed. We refer to that publication for all 
the details of the method and we only remind the reader of the basic formula, from 
which the RGI b-quark mass M^, can be implicitly derived: 

K-l 

LomB = LoTiLo, Mb) + 2-(''+^V^K) + LqAE , (4.15) 

k=0 

where the use of a finite volume scheme like the Schrodinger functional is assumed. 
In eq. (|4.15j) 

• Lq is the linear extent of the small volume where the matching between HQET 
and QCD is performed. 

• mB is the physical (spin-averaged) B^g^-meson mass. 

• r is an energy defined in terms of heavy-light correlators (with a heavy quark of 
mass M) in a finite volume L^. A relevant example is 

T{L,M) = i(rp + 3rv), Tp = ^ln[/A(xo-a)//A(xo + a)] , (4.16) 

{xq/T = 1/2 with T/L fixed) where Ty is defined analogously in the vector 
channel. Their static version rstat(-^) was defined in eq. (|3.15|) . Its infinite volume 
limit gives the static binding energy Egtat in eq. 1)2. 91) . 

• dm is a step scaling function 

a^{g\L)) = 2L [r,t.t{2L) - T.t^L)] . (4.17) 

Notice that in the difference the dependence of Fstat on 5m\\i cancels. Once the 
effective theory and the full one are matched in small volume the step scaling 
function provides, within the effective theory, the connection to large volumes, 
where contact with phenomenology can be made. This large scale difference is 
covered in several steps, such that at each stage scales differing only by a factor 
two appear. In this sense the sum in eq. 1)4. 15(1 connects the size Lq to = 2^ Lq, 
therefore Uk = g^{2^LQ) in the equation. 

• AE is an energy difference 

AS = ^,tat-rstat(ii^) . (4.18) 
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It is important to remark that in the way eq. ()4.15|) has been written ah the quantities 
on the r.h.s. are independent of dm-^ and can be computed in the continuum Hmit. 

In real applications of the method (see [9,36]) the binding energy -Estat has been 
computed on a lattice of size 1.6 fm while for Lk it sufficed to have almost half of 
it, more precisely Lk = L^ei = l-436ro. In view of combining this strategy with the 
use of the new static actions we have calculated the quantity aFgtat (do , -^^ref /«) > on the 
configurations generated in runs I- IV of table El The results are collected in table 01 





c-A 


„stat 
^A 

OS 


= 

^HYPl 


QHYP2 


6.0219 
6.1628 
6.2885 
6.4956 


0.5868(5) 
0.5575(4) 
0.5326(3) 
0.4948(3) 


0.5000(4) 
0.4763(4) 
0.4558(3) 
0.4243(3) 


0.2044(4) 
0.1921(4) 
0.1812(3) 
0.1643(3) 


0.1744(4) 
0.1625(4) 
0.1537(3) 
0.1390(3) 


P 


qA 


„stat 
^A — 

OS 


1-loop 

oHYPl 


QHYP2 


6.0219 
6.1628 
6.2885 
6.4956 


0.5871(5) 
0.5576(4) 
0.5335(3) 
0.4958(3) 


0.5003(4) 
0.4763(4) 
0.4567(3) 
0.4254(3) 


0.2050(4) 
0.1923(4) 
0.1822(3) 
0.1653(3) 


0.1769(4) 
0.1641(4) 
0.1547(3) 
0.1396(3) 



Table 3: Results for aT^ 



qo,Lrct/a) with L^d = 1.436ro, /3 = 6/po- 



The data can be described within one standard deviation by polynomial expressions. 
Explicitly, for the 1-loop value of cj^* and for 6 < /3 = Q/qq < 6.5 we provide the 
parameterizations {x = (3 — 6) 



arstat(s'o,iref/a) = 0.5916 

arstat(5o,ircf/a) = 0.5040 

aTsti,tigo,Lrc{/a) = 0.2068 

arstat(5o,iref/a) = 0.1787 



0.2140 X + 0.0418 j;' 
0.1726 X + 0.0285 x'^ 
0.0890 X + 0.0106 x'^ 
0.0915 X + 0.0255 x'^ 



for 

for SI 

for S^''' 
for 5?YP2 



(4.19) 
(4.20) 
(4.21) 
(4.22) 



with an uncertainty of about 2.5%o. In addition we see from the table that cj^* has very 
little impact on this quantity. 



5 Summary and conclusions 

We have proposed and investigated alternative discretizations for static quarks on the 
lattice. All of them have automatic 0(a) improvement, i.e. energies approach the con- 
tinuum limit with correction O(a^) if the light-quark sector is improved. The purpose of 
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introducing these actions was to achieve a better noise to signal ratio at large Euclidean 
times. This goal could be reached, since the exponential growth of the noise to signal 
ratio is non-universal and can be reduced by choosing a discretization with a small self 
energy for the static quark (relative to the lattice potential at the origin). The two 
actions with HYP links are best in this respect. 

Of course a prime criterion for choosing an action is to have small lattice artifacts. 
We have therefore investigated several quantities (figures EEEI)- It turns out that lattice 
artifacts are mostly indistinguishable for the different actions except for the case of 
Em, figure (and to a smaller extent ^a)) where the lattice artifacts with the two 
actions with HYP links are comfortably small compared to the ones with the standard 
action. Since an implementation of two or three static actions is no problem in any 
practical computation, it is advisable to compute with both S^^^'^ and S^^^^ and 
to perform continuum extrapolations constrained by universality. This promises to 
stabilize continuum extrapolations in a non-trivial manner, since e.g. figure El provides 
an example where the a-effects are different for the two actions. 

In view of future applications we have computed all coefficients needed to improve 
and renormalize the static-light axial current for all actions considered. Here, the renor- 
malization factors are now known non-perturbatively (in the quenched approximation) 
and the improvement coefficients are given up to uncertainties 0{gQ). In the course 
of these computations we again observed that the a-effects (now the ones before 0(a) 
improvement) are different for the different actions. 

We finally emphasize that in [43] it was shown that the force between static quarks, 
computed e.g. through Polyakov loop correlation functions, is automatically 0(a) im- 
proved. This is derived from the corresponding property of the Eichten Hill action for 
static quarks. This proof is also valid when the alternative actions are employed. Of 
course a significant gain in the noise to signal ratio can be achieved in this way as has 
first been seen in [29]. Such a gain in precision is especially relevant in the theory with 
dynamical quarks, where the noise-reduction methods of [44] are not applicable. In- 
deed, with the action S^, it has been possible very recently to observe string breaking 
in TVf = 2 QCD [45]. 

In addition, these discretizations can be used to efficiently calculate the 1/mb cor- 
rections to the static limit. In this context a preliminary study has been presented 
in [46]. 
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Table 4: The heavy-quark gluon vertex. 
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1 


-2if sin[a(p,(,) + ^)] 



Table 5: The 4 point heavy-quark gluon vertex. 



A Perturbation theory for the APE action 
A.l Feynman rules 

We use the conventions of [21,23]. Here we generalize the Feynman rules given for the 
EH action in appendix B.l of [21] (again we use the notation of that reference), for the 
action 

gA,a . w^'''{x,0) = {l-a)U{x,0)+aV{x,0) , (A.l) 

which reduces to eq. 1)2. 20() . when q = 1. We set 5mw = 0. The propagation of a 
static quark from the boundary of the SF to the bulk is given by the matrix iJ^ with a 
perturbative expansion 

oo 

H^{x) = ^g',Hl'\x). (A.2) 
fe=o 



20 



The terms proportional to with A; = 0, 1, 2 are given by 

= P+, (A.3) 

XO 6 

= -^EE^'^"E^J'^(P)«Mi)(^o-a + a.(z);p)r'^P+, (A.4) 

p uo=a 1=1 

hI'\x) = Hl'\xh + hI'\x)2 , (A.5) 

XQ UQ—a 6 
p.q uo=a vo=a i=0 

E 4'^(q)4o)(^o - a + a .(j); q)r'^r''P+ , (A.7) 

j=0 

XO 15 

J^f^(x)2 = ^j;X^e'(P+'i)-j;/if (p,q)g-(,)(no-a + a.(i);p) (A.8) 

p,q Mo=a j=0 

tii)i^o -a + a tii); q)T''T''P+ , (A.9) 
where the gluon field in the time momentum representation is defined by 

= :^5]e*P"go(xo,p), (A.IO) 
p 

^fcC^) = r^E«'^''"^"'"^'^9'^(^0'P) with A; = 1,2,3. (A.ll) 
p 

The vertex functions h^^^ and h^'^^ are listed in tables |1] and [51 They are Hnear in a and 
reduce to the ones for the EH action for a = 0. It is thus sufficient to give results for 
a = 1, to which we specialize from now on. 

A. 2 Self energy 

Here we compute, at one loop order, the linearly divergent contributions to the static 
propagator and the potential at zero distance for the action S^. We expand the corre- 
lation function fi^ix^) defined in eq. 1)3.13(1 

^/nx3)=i+5ovr^^^(x3,f) + .... (A.12) 

The Feynman diagrams that contribute to f^^'^^^ are given in fig. B.l of [35]. We then 
define 

As explained in sec. 2.1 the relevant quantity for the signal to noise ratio is Egeii — 
V{0)/2. In order to compute the divergent contribution of V{0) we introduce 

/j^^x) = -y ^ (Ch'(x)75Ch'(0) Ch(y)75Ch(z)) , (A.14) 
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which is related to /f^ via /i^'^(a;3) = {ajV)^ /|^'^(x). With the by now famihar 

notation for the one- loop coefficients, we then define 



lim f^^^'^^^ I X 
f-0 



^ a \ a 



The combination 



.(1) - 



(A.15) 



(A.16) 



is listed in tabled 

While the Feynman rules given above are valid for the result for r^^^ is the 
same for the action S\^. These two actions differ only by the normalization term 



f +Qtryt(^,0)y(x,0)) 



1/2- 



(A.17) 



multiplying the links W . This normalization factor is gauge invariant. Consequently, 
at one-loop order, the gluon fields which are obtained from the expansion of t, are 
connected only to themselves. There is no connected graph from t to the rest of the 
variables. In other words, these are tadpole graphs, which factorize (at one-loop!) in 
the sense that 



(Oh)t = (Oh)t=i X 



1 + 



(A.18) 



Here we have written down the expression for Oh which contains a single heavy quark, 
otherwise several terms oc t\ would appear. The function t\ is the result of the sum over 
the loop momentum of the tadpole graphs. Due to translation invariance in space it 
depends only on the time coordinate xq. The sum over xq in eq. ()A.18|) extends over all 
timeslices over which the heavy quark propagates. If, for the simplicity of the argument, 
we assume that we also have translation invariance in time, t\ does not depend on xq 
and is a pure (dimensionless) number (note that we do not have a factor a in front of the 
sum). It is then evident that e^"^^ is shifted by an amount ti, while SV^^^(^^ is shifted 
by 2ti and r^^^ is unchanged as claimed. 

As (at one-loop order and with periodic boundary conditions) the effect of t can 
entirely be absorbed into a change of e*^^\ it is also clear that the improvement coef- 
ficients c^*^* , 6^*^* are independent of t and in particular they are identical for and 
SI. 



A. 3 Improvement coefficients 

To compute the improvement coefficients of the static-light axial current we make a 
one loop perturbative expansion of a correlator involving the static-light 0(a) improved 
axial current 



/IT (^o) = !T\^^) + a cr/ir (^0 



(A.19) 
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The perturbative expansion reads 

fri^o)=^9i'fT'^'\xo), (A.20) 

A:=0 



and the Feynman diagrams contributing to the one loop term are summarized 

in figure B.l of [21]. 

In the following, if not explicitly indicated, we will insert the local operators always 
at xq = T/2 and supress that argument. 

A. 4 Determination of c^^* 

The universality of the continuum limit gives us a handle to compute the improvement 
coefficient for a certain static action, once the improvement coefficient is known for 
another static action. This remark, given two actions with the same continuum limit, 
is valid in general and not restricted to the static case. It can be translated into the 
fact that the ratio 

ftfiSh) 



Pi{s^;e,e')= 'Z , (A.21) 

Ja,i Wh)\e' 

is independent of the action used up to O(a^), once suitable values of the improvement 
coefficients are chosen. The basic idea is then to use the already known value of 
for the EH action [21,34] and to enforce 

pjiS^''-e,9')=piiSt,0,O') , (A.22) 



in order to compute (f^^^'^^^ for S^. Different choices for 9 and 9' will give equivalent 
definitions for the improvement coefficient, differing only by cutoff effects 0{a/L). The 
perturbative expansion of eq. (|A.21() reads 

oo 

Piis^;e,9') = Y,9fpf\s^-.e,9') . (A.23) 

fe=0 

After some trivial algebra we obtain 

pj') (5h; e, 9') = (5h; 9, 9') + ac^'^^^ (S^) [a{9) - A{9')] , (A.24) 

and 

j.stat,(l) / r, XI j,stat,(l) , Ni 

P"{^h,0,9}- — — , (A.25) 

Ja \e Ja \9' 

where we use the fact that the correlator at tree-level does not depend on the static 
action, and 

nstat,(0) I 

Ja \d 
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An explicit formula for the desired improvement coefficient is 

stat,(l)/„A\ _ stat,(l)/crEH\ 



>{S^) = cr^'-'iSn- lim Ac 

(A.27) 



Ac 



pW(st,e,e')-pWis^^;e,e') 

A{9) - A{e') 



We have performed a one loop computation of Ac for 9, 9' € {0.0, 0.5, 1.0}, T = L, and 
for resolutions L/a = 4, 6, ...,48. The bare light quark mass was set to the critical 
mass 

= m«g2 + Oig^) , (A.28) 

with [47] 

amW = -0.2025565 x Cp . (A.29) 

The input value 

^^tat,(l)^^EH) = -i- X 1.0351(1) , (A.30) 
was taken from [34]. In figure|Slwe plot Ac as a function of a/L for the 3 possible choices 



-0.07 
-0.08 
-0.09 
-0.1 
-0.11 

0.1 0.2 (a/L) 

Figure 8: Dependence of Ac, eq. fA.27|) . on a/L. 

of {9,9'). It is clear that the limit a/L — > can be taken with reasonable precision. The 
plot indicates also that defining the improvement coefficient without performing the 
continuum limit will give a perturbative 0(a/L) uncertainty of the order of 0.02 g^. 
To obtain the final number eq. (|3.4jl we used both the extrapolation procedures of [48] 
and [49], finding consistent results. 
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A. 5 Determination of U 



A 



To determine 6^*^* we follow the strategy just applied to compute c^^*. The ratio 



pi{Sh;m' 



(1 + abr (^h)S)/A.f(>gh)|> 
{l + abT\Si,)m^)r^f{St)[, 



(A.31) 
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has a well defined continuum limit. Once suitable values of the improvement coefficients 
are chosen, the continuum limit is approached with a rate of O(o^), as long as one 
considers (i) the quenched approximation or (ii) is interested in one-loop accuracy, only. 
These restrictions are necessary, since for full QCD and starting at order Qq, additional 
terms proportional to aniq with coefficients denoted by b(^^^,bg (see [21]) are needed to 
cancel all 0(a) effects (the term proportional to 6^ cancels trivially in the above ratio). 
We can hence require 



~ / cEH / 



pi(5h ;mq,mq) + 0(5, 



0) ' 



which yields immediately 



, stat,(l) 



Ab 



,stat,(l)^C'EH 
"A 



(5^)+ lim Ab, 

a/L-,0 



(A.32) 



(A.33) 



a(m^ 

/.stat,(l) 

/a,i 



/.stat,(l) 
JA,l 



i.stat,(l) 
^A,I 



{Sh 



EH^ 



j,stat,(0) 
J A 



j,stat,(0) 
J A 



^stat,(0)(^A)|^ 



and choose m'^ such that 



In the numerical evaluation, we set 

Lm'^ = L(Zmmq(l + bmam'^)) = 0.24 . 



(A.34) 



(A.35) 



Here, in contrast to the MC evaluation, has been defined in the lattice minimal 
subtraction scheme, where = 1 + ^ln{L/a)gQ. The improvement constant 6^ = 

— ^ — 0.07217 Cf Qq is known from [50] and the input value of 5^^*'*-^^ for the EH action 



-0.056(7) , 



(A.36) 



was computed in [21]. 

In figure IHl we plot Ab as a function oi a/L for the 3 possible choices of 6. Again 
the 0{a/L) effects are well under control and almost independent of the choice of 9. 
Extrapolating a/L ^ as above yields the result eq. (|3.8j) . 



B Approximate one-link integral 

The basic idea of the one-link integral is that (in the pure gauge theory) correlation 
functions of Polyakov loops remain unchanged if the link variable U (x, 0) is replaced by 



U{x,0) 



fdUU exp {12/g^ Re tr [UV^x, 0)]) 



(B.l) 



/ dC/ exp (12/^2 Re tr[UV^x,0)]) ' 

whereas the variance (and therefore the error) is substantially reduced. While U{x,0) 
would thus be a possible choice for W{x,0), its expression is rather unhandy for gauge 
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Figure 9: Dependence of Ab, eq. fA..S4|) . on a/L. 



(a/L) 



groups SU(A^) with N > 2, [51,52]. A more convenient approximation can be inferred 
form the SU(2) case, where the group integral is proportional to V{x,0) and takes the 
form 

Vl2{24v/gi) 



(B.2) 



in terms of the Bessel functions Ii and l2- For SU(3) we make a similar ansatz 

W{x,0) = g^^V{x,0)sig^''v) . (B.3) 

Although this can not reproduce the SU(3) one-link integral exactly, since that is a 
more complicated function of V and V"^, it turns out that the simple form 



s(y) = fc/(y + l/5) , 



(B.4) 



with k = 1.092 comes very close and thus is a good candidate to define an action with 
reduced noise and roughly unchanged signal with respect to the Eichten-Hill action. 
The form of s{y) has been obtained numerically by minimizing the difference between 
the norm of the links U{x,0) (constructed by a multihit procedure) and the norm of 
the links in eq. (|B.3() with s{y) replaced by some trial functions of y. For this tuning 
we used a set of quenched configurations with SF boundary conditions on an 8^ lattice 
at /3 = 6 and /? = 6.5. ^ 



C Simulation parameters 

We list in table IHl the lattice sizes and the bare parameters of our simulations. The 
angle defines the periodicity of the fermionic fields in all the spatial directions. We set 
T/L = 1 and choose the background field to vanish. The first n value always corresponds 

^The actual value of k is irrelevant for physics as it only influences the self energy. We chose k = 1 
in eq. (EHJ- 
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run 


L/a 


P 


K 





7f coni. 


I 


8 


6.0219 


0.135081, 0.1344011 


0.0, 0.5, 1.0 


4800 


II 


10 


6.1628 


0.135647, 0.1351239 


0.0, 0.5, 1.0 


3520 


III 


12 


6.2885 


0.135750, 0.1353237 


0.0, 0.5, 1.0 


4000 


IV 


16 


6.4956 


0.135593, 0.1352809 


0.0, 0.5, 1.0 


3200 


ZI 


6 


6.2204 


0.135470 


0.5 


4000 


sZI 


12 


6.2204 


0.135470 


0.5 


7200 


ZII 


8 


6.4527 


0.135543 


0.5 


3200 


sZII 


16 


6.4527 


0.135543 


0.5 


3200 


ZIII 


12 


6.7750 


0.135121 


0.5 


7360 


sZIII 


24 


6.7750 


0.135121 


0.5 


2000 



Table 6: Collection of the simulation parameters and statistics. 



to Kc, where the PCAC mass vanishes [35,42], while the second one for the runs I to 
IV has been chosen such that Lmn = 0.24 with mR, the running quark mass in the SF 
scheme at the scale = 2Lmax = 1.436ro: 

mB. = Zmfhq, mq = mq(l + 6ma"^q) , «"iq=^(-- — ] . (C.l) 

The low energy scale rg (— 0.5 fm) has been introduced in [37] and in [38] the ratio rg/a 
has been computed to a precision of about 0.5% in a range of /3 values which includes 
runs I to IV in table |HI For the coefficients and 6jn in eq. (|(].1|1 we have used the 
results in [39] and the improvement coefficients Ct,Ct were set to their two-loop and 
one-loop values respectively [19,40,53]. 
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